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ABSTRACT 

Transverse oscillations of solar filament and prominence threads have been fre- 
quently reported. These oscillations have the common features of being of short pe- 
riod (2-10 min) and being damped after a few periods. Kink magnetohydrodynamic 
(MHD) wave modes have been proposed as responsible for the observed oscillations, 
whereas resonant absorption in the Alfven continuum and ion-neutral collisions are the 
best candidates to be the damping mechanisms. Here, we study both analytically and 
numerically the time damping of kink MHD waves in a cylindrical, partially ionized 
filament thread embedded in a coronal environment. The thread model is composed of 
a straight and thin, homogeneous filament plasma, with a transverse inhomogeneous 
transitional layer where the plasma physical properties vary continuously from fila- 
ment to coronal conditions. The magnetic field is homogeneous and parallel to the 
thread axis. We find that the kink mode is efliciently damped by resonant absorption 
for typical wavelengths of filament oscillations, the damping times being compatible 
with the observations. Partial ionization does not affect the process of resonant absorp- 
tion, and the filament plasma ionization degree is only important for the damping for 
wavelengths much shorter than those observed. To our knowledge, this is the first time 
that the phenomenon of resonant absorption is studied in a partially ionized plasma. 

Subject headings: Sun: oscillations — Sun: magnetic fields — Sun: corona — Sun: 
prominences 



1. INTRODUCTION 



The fine- structures of solar prominences and filaments are clearly seen in high-resolution ob 
servations. These fine- structures, here called threads, appear as very long (5" - 20'^ ) and thin 
(0".2 - 0".6) dark ribbons in Ha images of filaments on the solar disk (e.g.. lLinll2004l : iLin et al. 
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2005L 120071 12008L 120091) . as well as in observations of promi nences in the solar limb from the 



Solar Optical Telescope (SOT) aboard the Hinode satellite (e.g. JOkamoto et al.ll2007l: iBerger et al 



20081 : IChae et alJUOOSl : iNing et all 120091) . Although statistical studie s show that the orientation 
of threads can significantly vary within the sa me filament jLinI 120041). vertical t hreads are more 
commonly seen in quiescent prominences (e.g., Berger et al. 20081 : Chae et al. 2008i) whereas hori- 



zontal threads are usually observed in active region prominences (e.g.Jo kamoto et aPboOV ). From 



the theoretical point of view, fil ament threads have been modeled as ma gnetic flux tubes anchored 
in the solar photosphere (e.g., iBallester & Priest! 1 19891 : iRempel et al.l ll999). In this interpreta- 
tion, only part of the flux tube would be filled with the cool (~ 10"^ K) filament material, which 
would correspond to the observed threads. It has been also suggested by differential emission mea- 
sure studies that each thread might be surrounded by its own prominence-corona transition region 
(PCTR) where the plasm a physical properties would abruptly vary from filament to coronal con- 
ditions (^ Cirigliano et al. 2004) . The filament material, roughly composed by 90% hydrogen and 
10% helium, is only partially ionized for typical filament temperatures, although the precise value 
of the ionization degree is not well-known a nd could probablv varv in different filaments or even 
in different threads within the same filament (IPatsourakos & Vialll2002h . 



Oscillations of prominence and filament threads have been frequently rep orted since tele- 
sco pes with a high time and spatial resolution became available. Early works by lYi et al.l (1199 1|) 
andiYi & Engvoldl (|l99ll), with a relatively low spatial resolution (~ 1"), detected oscillatory vari- 
ations in Doppler signals and He I intensity from threads in quiescent filaments. Later, Ha and 
Doppler observations with a much better spatial resolutio n (~ 0.'^ 2) found evidence of osci llations 
and propagating waves along quiescent filament threads (lLinll2004l : iLin et al.l 120071 12009|) . while 
observations from the Hinode spacecraft showed transv erse oscillations o f threadlike structures in 
both active region dOkamoto et al.ll2007b and quiescent (INing et al.ll2009l) prominences. Common 
features of these observations are that the reported periods are usually in a narrow range between 
2 and 10 minutes, and that they are of small amplitude, with the velocity amplitudes smaller than 
~ 3 km s~' . Some theoretical works have attempted to explain these observed oscillations in terms 
of linear magnetohyd rodynamic (MHD) waves supported by the th read body, modeled as a cylin- 
drical magnetic tube (|Diaz et al.ll2002l : iDymova & Rudermanll2005[) . These studies concluded that 
the so-called kink MHD mode is the best candidate to explain transverse, nonaxisy metric thread 
oscillations, and is als o consistent with the r eported short periods. An applic ation of this interpreta- 
tion w as performed bv lTerradas et al.l(l2008b. who made use of the model by lOymova & Ruderman 
and the observati ons by [O kamot o et al.l (120071) to obtain lower limits of the prominence 
Alfven speed. Similarly. ILin et al.i (2009) interpreted their observations of swaying threads in Ha 
sequencies as propagating ki nk waves and gave a n esti m ation of the Alfyen speed. The reader i s 
refereed t o rece nt reviews by lOliver & Ballesteil (|2002l) : iBallesteil (|2006|) : iBanerjee et al.l (|2007|) : 
Engvoldl (|2008[) for more extensive comments. 



-3- 



Another interesting characteristic of prominence oscillations is that they seem to be damped 



after a few periods. A : 



Landman et al 



Molownv-Horas et al. 



though this behavior was previou sly suggested by the results of some works 



1977 



Tsubaki & Takeuchi 



1986 ). it was first extensively investigated by 



(|1999h and iTerradas et al.l (|2002|) . These authors studied two-dimensional 



Doppler time-series from a quiescent prominence and found that oscillations detected in large ar- 
eas of the prominence were typically damped after 2-3 periods. Similar results were obtained in 
a more recent work by lMashnich et al.l (120091) . This damping patter n is also se en in several high- 
resolution Doppler time- serie s from individual filament threads by iLinI (120041) . as well as in the 
Hinode/SOT observations by iNing et al.l (|2009h . who reported a maximum number of 8 periods 
before the oscillations disappeared. In the context of the kink MHD n iode interpret at ion, several 



mechanisms have been proposed to explain this quick damping (see IOlivenl2009|) . ISoler et al 



(120081) studied the damping by nonadiabatic effects (radiative losses and thermal conduction) in a 
fully ionized, homogeneous, cylindrical filament thread embedded in a homogeneous corona, and 
obtained a kink mode damping time larger than 10^ periods for typical prominence conditions, 
meaning than nonadiabatic effects cannot explain the observed damping of transverse oscillations. 
Subsequently, lArregui et al.l (l2008b considered a transverse inhomogeneous transitional layer be- 
tween the fully ionized filament thread and the corona, and investigated the kink mode damping 
by resonant absorption in the Alfven continuum. They obtained a dampi ng time of approxim ately 
3 periods for typical wavelengths of prom inence oscilla t ions. L ater on, ISoler et al.l (|2009al . here- 
after Paper I) complemented the work by lArregui et al.l (|2008r) by also considering the damping 
by resonant absorption in the slo w continuum, but con cluded t hat this effect is no t relevant and 
obtained similar results to those of lArregui et al.l (|2008|) . Finally. ISoler et al.l (|2009bl hereafter Pa- 
per II) included the effect of partial ionization in a homogeneous filament thread model. These 
authors found that ion-neutral collisions can efficiently damp the kink mode for short wavelengths, 
although they are not efficient enough in the range of typically observed wavelengths. 

On the basis of these previous studies, resonant absorption seems to be the most efficient 
damping mechanism for the kink mode and the only one that can produce the observed damping 
times. On the other hand, the effect of partial ionization could be also relevant, at least for short 
wavelengths. The aim of the present work is to take both mechanisms into account and to assess 
their combined effect on the kink mode damping. To our knowledge, this is the first time that the 
phenomenon of resonant absorption is studied in a partially ionized plasma. The filament thread 
model assumed here is similar to that of Paper I and is composed of a homogeneous and straight 
cylinder with prominence-like conditions embedded in an unbounded corona, with a transverse 
inhomogeneous transitional layer between both media. We consider the prominence material to be 
partially ionized. As in Paper II, we u s e the l inearized MHD equations for a partially ionized, one- 
fluid plasma derived byl Forteza et al.l (|2007|) . Since we focus our investigation on the kink mode, 
we adopt the/? = approximation, where p is the ratio of the gas pressure to the magnetic pressure. 
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This approximation removes longitudinal, slowlike modes but transverse modes remain correctly 
described. In the next Sections, we investigate the kink mode damping by means of both analytical 
and numerical computations, and assess the efficiency of the considered damping mechanisms. 

This paper is organized as follows: Section [2] contains a description of the model configura- 
tion and the basic equations. The results are presented and discussed in Section |3l Finally, our 
conclusions are given in Section IH 



2. MODEL AND METHOD 



2.1. Equilibrium Properties 



We model a filament thread as an infinite and straight cylindrical magnetic flux tube of radius 
a surrounded by an unbounded coronal environment. Cylindrical coordinates are used, namely 
r, (f, and z, for the radial, azimuthal, and longitudinal directions, respectively. In the following 
expressions, a subscript indicates equilibrium quantities while we use subscripts / and c to 
explicitly denote filament and coronal quantities, respectively. The magnetic field is homogeneous 
and orientated along the cylinder axis, = B^e^, with Bq = 5 G constant everywhere. We 
adopt the one-fluid approximation and consider a hydrogen plasma composed of ions (protons), 
electrons, and neutral atoms. Since we are interested in kink MHD waves supported by the filament 
thread body, we also assume the /3 = approximation, which neglects gas pressure efi'ects and 
removes longitudinal, slowlike modes. With such assumptions (see details in lForteza et al.l 120071 
and Paper II), the plasma properties are characterized by two quantities: the fluid density, po, and 
the ionization fraction, JIq, which gives us information about the plasma degree of ionization. The 
allowed values of JIq range between fxo = 0.5 for a fully ionized plasma and //o = 1 for a neutral 
plasma. 

The density profile assume d here only depends on the r adial direction and is the same as in 
Paper I (which was adopted after IRuderman & Robertsll2002h . namely 



with 



Po {r) = 




if a- 1/2 < 



r<a- 111, 

r<a + l/2, 
r>a + 1/2, 



1 + — - 1 sm 



(1) 



(2) 



Pf/ \ Pfi 

where we take p/ = 5 x 10"" kg m"^ and pc = 2.5 x 10"'^ kg m"^, the density contrast between 
the filament and coronal plasma being Pf/pc = 200. In these expressions, a is the thread mean 
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radius and I is the transitional layer width. The limit cases Ija = and l/a = 2 correspond to 
a homogeneous thread without transitional layer and a fully inhomogeneous thread, respectively. 
We also take the same functional dependence for the ionization fraction profile, 



p-f, if r < a - 1/2, 

P-o W = { P-tr (r) , if a - 1/2 < r < a + 1/2, 
pc, if r > a + 1/2, 



(3) 



with 



j(r-a) 



(4) 



where the filament ionization fraction, p f, is considered a free parameter and the corona is assumed 
to be fully ionized, so pc = 0.5. 



2.2. Basic Equations 



We co nsider the basi c MHD equations for a partially ion ized plasma in the one-fluid approach 
derived by iForteza et al.l (|2007|) (see also iBraginskiil Il965h . Since we are dealing with small- 
amplitude perturbations, we restrict ourselves to linear perturbations. After removing gas pressure 
terms by setting /3 = 0, the relevant equations for our investigation are the momentum and the 
induction equations, whose linearized version in MKS units are. 



Po^ = -[(VxBi)xBo], 
ot n 



(5) 



^Bi 

~5r 



= Vx(vi X Bo)-Vx(77V X Bi) + Vx{77a [(V x Bj) x Bq] x Bo}-Vx[77h (V x BJ x Bo] , (6) 



along with the condition V • Bi =0. In these equations Bi = {Br,B^,B^, Vi = {yr,v^,v^ are 
the components of the magnetic field and velocity perturbations, respectively, while = 4;r x 
10"^ N is the vacuum magnetic permeability. Quantities rj, r]A, and t/h in Equation Q are the 
coefficients of ohmic, ambipolar, and Hall's magnetic diffusion. Note that these nonideal terms of 
the induction equation appear due to collisions between the different plasma species. In particular, 
ohmic diffusion is mainly governed by electron-ion collisions and ambipolar diffusion is mostly 
caused by ion-neutral collisions. Hall's effect is enhanced by ion-neutral collisions since they tend 
to de couple ions from the ma gnetic field while electrons remain able to drift with the magnetic 
field (|Pandey & Wardlell2008|) . The ambipolar diffusivity is commonly expressed in terms of the 
Cowling's coefficient, rjc, as follows, 

= ^0 ■ (V) 
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As explained in Paper II, rj and ric correspond to the magnetic diffusivities longitudinal and perpen- 
dicular to magnetic field lines, respectively. Due to the presence of neutrals, rjc » ?7, meaning that 
perpendicular magnetic diffusion is much more eflicient than longitudinal magnetic diffusion in a 
partially ionized plasma. In a fully ionized plasma rjc = rj, so magnetic diffusion is then isotropic . 
Since rj, nc, and 77^ are functions of the pl a sma physical conditions (see, e. g., 'Braginskii 1965 : 
2OO4I : Leake et allllOOsI : Leake & Arbeilbood IPandev & Wardle 20081 Paper 



Khodachenko et al. 



II for expressions), their values in our equilibrium depend on the radial direction. It is convenient 
for our following analysis to write ohmic. Cowling's, and Hall's diffusivities in a dimensionless 
form. 



T] = 



m = 



(8) 



VA/a VA/a VA/a 

where the tilde denotes the dimensionless quantity and va/ = BqI ^n^Pf is the filament Alfven 
speed. Note that ff^ is usually called the magnetic Reynolds number. Figure [U displays the ra- 
dial profiles of f), fjc, and t/h according to the equilibrium properties. The Figure focuses in the 
transitional layer, where the dimensionless coefficients vary by several orders of magnitude from 
internal to external values. 
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Fig. 1. — Radial profiles of (a) 77, {b) fjc, and (c) ?7h considered in this work. The transitional 
layer (shaded zone) is enhanced in order to see the change from filament to coronal values. The 
linestyles represent different ionization degrees: fxf = 0.5 (dotted line), Mf = 0.6 (dashed line), 
fif = 0.8 (solid line), and /i/ = 0.95 (dash-dotted line). 

Since (p and z are ignorable coordinates, perturbations are put proportional to exp (imcp + ik^z - icot), 
where oj is the oscillatory frequency, and m and are the azimuthal and longitudinal wavenumbers, 
respectively. Then, Equations ([S])-® become, 

coVr = -'^(KB, + iB'y (9) 
Co 

v\ /m \ 
ajv^ = j-[-B,-LB^j, (10) 
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ajv, = 0, (11) 
coBr = -BoKvr + r]\-B'^ + -B^ - i—BA + t/c (KB[ - iklBr) + 7/h5o [iK-B, - ik^B^j , (12) 

loB^ = -BoLv^ + 7] + i^B'^ - i^B^ + ^B'^ - + 77' + i^B^ + 

+77C {ih'^B, - iklB^j + 77h5o [ik^Br - KB[) , (13) 

(\ in \ I 1 wi^ 1 tn \ 

iv' + i-Vr + '/c ^'^r + 2-5- - r5j. + + k^-B,- + ik^—B^ 
r r j y i y - yi y r j 

+r]'c (KBr + iB'^ + 77h5o + k}-B^ - + r]'^Bo (k,B^ - ^5,) , (14) 

where the prime denotes derivative with respect to r. From Equation (fTT)) we get v,- = 0, so no 
longitudinal displacements are allowed in the yS = approximation. These equations form an 
eigenvalue problem, with co the eigenvalue and (v^, v^, B^, B^, B,^ the eigenvector. In the present 
work, we study the time damping of kink waves, so we assume m = 1 and take a real and positive 
A complex frequency, co = cor + icoi, is therefore expected, damped solutions corresponding 
to oji < 0. The period, P, and damping time, r^, are related to the real and imaginary parts of the 
frequency as follows, 

P = -. TO = ^. (15) 



2.2.1. Importance of Hall's Term 

In Paper II, we showed that ohmic diffusion can be neglected in front of Cowling's diffusion 
for large kya, while ohmic diffusion dominates for small values of k^a. In the present work, we 
have also included Hall's diffusion. Here, we compare the importance of Hall's term with the other 
nonideal terms in the induction equation. To achieve this, we first have to assess the typical length- 
scales of the dissipative terms of the induction equation. We note that longitudinal derivatives 
only appear in the terms with fjc and rju of Equations (fT2l) -(fT4l). while the terms with fj only 
contain radial and azimuthal derivatives. Therefore, the typical length-scale of both Cowling's and 
Hall's terms is the longitudinal wavelength, A^, which can be expressed in terms of the longitudinal 
wavenumber, = 2nk~^. On the other hand, the typical length- scale of the ohmic term is the 
filament thread mean radius, a. By taking into account these relevant length- scales, we define the 
dimensionless number 'H as the magnitude of Hall's term with respect to that of the ohmic term, 

|Vx[;7h(VxBOxBo]| (KA'm 

|Vx(?7VxBi)| ^\2;rj 77' ^ 
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From Equation (fT6l) . one can compute the value of k^a for which Hall's term becomes more impor- 
tant than the Ohm's term by setting 'H ~ \. So, one obtains, 




{k^a)^^ln ^\ — . (17) 



We can compare this transitional wavenumber with that obtained in Section 2.3. of Paper n that 
determines when Cowling's diffusion begins to dominate over ohmic diffusion. 




{k^a)r^2n.\^. (18) 



Since in our equilibrium fic > fju, one gets (k^a)^ < (k^a)^, meaning than Cowling's diffusion 
becomes dominant over ohmic diffusion for smaller wavenumbers than those needed for Hall's 
term to become relevant. Next, we have to know whether Cowling's diffusion or Hall's diffusion 
dominate for kM beyond both transitional values. To do so, we define the dimensionless number 
"Kc as the magnitude of Hall's term with respect to that of the ambipolar term, 

|Vx[?7h(VxBi)xBo]| ^ 
|Vx{^A[(VxBi)xBo]xBo}r ^c' 

Considering again that rjc > Uh, the condition "TYc < 1 is always satisfied, therefore Cowling's 
diffusion is always more important than Hall's diffusion. By means of this simple dimensional 
analysis, we expect the presence of Hall's term to have a minor effect on the results. 



2.3. Analytical Dispersion Relation 

Some analytical progress can be performed when no transitional layer is present, i.e., l/a = 0. 
We showed in Paper II that it is possible to give an analytical dispersion relation when the terms 
with rj and rju are dropped from the induction equation. In such a situation, the induction equation 
(Equation Q) can be written in a compact form as follows. 



— i = ^Vx(viXBo). 
dt vi 



(20) 



where = 



- iojric is the modi fied Alfven speed squared (iForteza et al.ll2008l). Next, we follow 



the standard procedure (see, e.g., lEdwin & RobertsI Il983l : iGoossens et al.l 120091) and obtain the 
dispersion relation of trapped waves by imposing that the radial displacement, = ivr/co, and 
the total pressure perturbation, = BoB^/fi, are both continuous at the thread boundary, r = a. 
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and that perturbations must vanish at infinity. So, the dispersion relation governing transverse 
oscillations is D,„ (o), = 0, with 



Dm (oj, k,) = 



He 



K' {n,a) 



nif 



(21) 



where 7^ and K,n are the Bessel function and the modified Bessel function of the first kind (lAbramowitz & Stegun 
1972t) . respectively, and the quantities nif and n,. are given by 



m. 



r2 



^ Ac 



(22) 



Note that Equation (121]) is valid for any value of m. 



Next, our aim is to extend this analytical analysis to the case l/a 0. When an inhomo- 
geneous transitional layer is present in the equilibrium, the kink mode is resonantly coupled to 
Alfven continuum modes. The radial position where the Alfven resonance takes place, ta, can be 
computed by setting the kink mode frequency, oj^, equal to the local Alfven frequency, ojp, = Vpjc^. 
The expression of ta corresponding to our density profile was obtained in Paper I (Equation (9)£|, 



I 

a -\ — arcsm 
n 



Pf+Pc 
Pf-Pc 



^^A/'^z 



Pf 



(pf-pc) 



(23) 



The ideal MHD equations are singular at r = r^. This singularity is removed if dissipative effects, 
such as magnetic diffusion or viscosity, are considered in a region around the resonance point, i.e., 
the dissipative layer. A method to obtain an analytical dispersion relation in the presence of an 
inhomogeneous transitional layer is to combine the jump conditions at the resonance point wit h 
the so-called thin boundary (TB) approxim ation, which w a s first u sed b vlHoUweg & Yand(ll988b . 
The jump conditions were first derived by ISakurai et al.| (Il991al) and iGoossens et al.l (1 1995b for 



the driven problem, and later by lTirry & GoossensI (|l996l) for the eigenvalue problem. They have 



been used in a number of papers in the context of MHD waves in the solar atrnosphere (e.g.. 



Sakurai et al 



1 99 ibi: IGoossens et al.lll992l: iKeppens et al.lll994l : iStenuit et al.lll998l : lAndries et al 



2OO0I : I Van Doorsselaere et al.l 



20041 . P aper I among o t her vy orks). An important result for the 



present investig ation was obtained b y IGoossens et al.l (Il995|) . who proved that the jump condi- 
tions derived by lSakurai et al.l (|l991ar) in i deal MHD rernain vali d in dissipative MHD. This allows 
us to apply the jump conditions derived by lSakurai et al.l (|l991ar ) to our case. Th e assumptions be - 
hind the TB approximation and its applications have been recently reviewed by IGoossens (2008). 



'Note that there is a typographical error in Equation (9) of Paper I since the term (which corresponds to ourp^) 
should be -pg. 
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in short, the main assumption of the TB approximation is that there is a region around the dissipa- 
tive layer where both ideal and dissipative MHD applies. In our equilibrium, we can assume that 
the thickness of the dissipative layer, namely 6p„ roughly coincides with the width of the inhomo- 
geneous transitional layer. This condition is approximately verified for thin layers, i.e., Ija 1. 
So, we can simply connect analitycally the perturbations from the homogeneous part of the tube to 
those of the external medium by means of the jump conditions and avoid the numerical integration 
of the dissipative MHD equations across the inhomogeneous transitional layer. 

The jump condition s for the radial displacement and the total pressure perturbation provided 



by lSakurai et al.l (Il991a|) in the case of a straight magnetic field are. 



m] = -i7:--^pj, [bT]]=0, at r = rA, (24) 

where [[X]] = Xc -Xf stands for the jump of the quantity X, and A = ^ (^a>^ - oj^^. By considering 
that in our model the magnetic field is straight and constant so that the variations of the local 
Alfven frequency are only due to the variation of the equilibrium density, we can write |poA|^^ = 
a>\\drPo\rj^. Ill addition, from the resonance condition we have a)\ = o*^. Applying the jump 



conditions (|24)) . we arrive at the dispersion relation in the TB approximation, 

2 / 2 

m /r. 

(CO, k,) = -in , ^ ^, , (25) 

with Dm {co, ky) defined in Equation (|2T]) . Note that in order to solve Equation (|25l) we need the 
value of the kink frequency, co/,. For this reason, we use a two-step procedure. First, we solve the 
dispersion relation for the case l/a = (Equation (|2T]) ) and obtain cok. Next, we assume that the 
real part of the frequency is approximately the same when the inhomogeneous transitional layer is 
included, allowing us to determine from Equation (12^ and therefore |<9^pol^^. Finally, we solve 
the complete dispersion relation (Equation (l25l) ) with these parameters. 



2.3.1. Expressions in the Thin Tube Limit 



Equation (1251) is a transcendental equation that has to be solved numerically. As in Paper I, it is 
possible to go further analytically by considering the thin tube (TT) approximation, i.e., k^a «; 1. 
The combination of both the TB and TT approximations has been done in several works (e.g., 
Goossens et aDll992l : lRuderman & Robertslbooilooossens et al.ll2002[l2009h . In the TT limit, we 
approximate the kink mode frequency as 



1 +Pclpf 



VhfK- 



(26) 
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We put this expression in Equation (l23l) and obtain that ~ a. Next we perform a first order, 
asymptotic expansion for small arguments of the Bessel functions of Equation (|25] ). The dispersion 
relation then becomes, 



I 2 ,2t2 \ I 2 ,2^2 \ . (m\ Pfpc (^^ - K^lf)(^^ - ^z^lc) 



(27) 



We now write the frequency as co = cot + and the modified Alfven speed squared is approxi- 



mated by r 



• icokrjc- We insert these expressions in Equation (ITTl) and neglect terms with a»j 



and ojik:. It is straight-forward to obtain an expression for the ratio coi/cok. 



n /m\ PfP, 
2U) 



(pf + p,) \drPo\ 
1 [pf7]Cf+pcr]Cc)K 
^ ^l(pf+Pc)(pfAf+PcA)j 



1 [Pf-Pc) 

4 Pfpc 



( 



1 



1 



— + - 

\<f <c> 



ricf-ricc 



(28) 



The first term in Equation (|28l) owes its existence to the factor in the dispersion relation related 
to the TB approximation and represents the contribution of resonant absorption. For long wave- 
lengths, the factor (l /v^^ -I- 1 lv\^ rjCfTlCc^ll'^ can be neglected, so the term related to the resonant 
damping is independent of the value of Cowling's difFusivity and, therefore, of the ionization de- 
gree. Then, one can see that in the TT case the term of Equation due to the resonant damp- 
ing takes the s a me fo rm as in a fully ion ized plasma, which was prev iously obtained by, e.g.. 



Goossens et al. 



1991 Equation (77)) and iRuderman & RobertsI (12001 Equation (56)). On the 



other hand, the second term in Equation (|28] ) is related to the damping by Cowling's diffusion and 
is also present in the case l/a = 0. This term is proportional to h, so we also expect it to be of a 
minor influence in the TT regime. 

Next, we take into account that for the present sinusoidal density profile, l^rPol^^ ~ n(^pf - pc) /2/. 
We insert this expression in Equation (|28] ) and then use it to give a relation for the ratio of the 
damping time to the period. 



P 



n 



I 

ml - 

a 



Pf-Pc 



2(p/77c/ +Pcfjcc)k,a 



(29) 



where we have used Equation ^ to rewrite Cowling's diffusivities in their dimensionless form. To 
perform a simple application, we compute td/P from Equation (|29l ) in the case m = 1, ha = 10"^, 
and l/a = 0.2, resulting in td/P « 3.18 for a fully ionized tread (Jif = 0.5), and r^/P « 3.16 
for an almost neutral thread (fif = 0.95). We note that the obtained damping times are consistent 
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with the observations. Furthermore, the ratio t^/P depends very slightly on the ionization degree, 
suggesting that resonant absorption dominates over Cowling's diffusion. To check this last state- 
ment, we compute the ratio of the two terms on the right-hand side of Equation (|29l ). which allows 
us to compare the damping times exclusively due to resonant absorption, (td)^^, and Cowling's 
diffusion, (td)c5 



(td) 



RA 



2(pf+ p,) / p^fjcf + Pcflcc \ Ka ^^^^ 
(td)c \| Pf \ Pf-Pc I m(l/ay 
This last expression can be further simplified by considering that in filament threads pf pc and 
fjcf^ fjcc, so that 

« V2^c,-^. (31) 
(td)c mil /a) 

Thus, we see that the efficiency of Cowling's diffusion with respect to that of resonant absorption 
increases with k^a and juf (through ^c/)- Considering the same parameters as before, one obtains 
(td)ra / (td)c ~ 2 X 10"^ for fif = 0.5, and (td)ra / (td)c ~ 6 x 10*^ for/// = 0.95, meaning that 
resonant absorption is much more efficient than Cowling's diffusion. From Equation (I3TI ) it is also 
possible to give an estimation of the wavenumber for which Cowling's diffusion becomes more 
important than resonant absorption by setting (td)ra / (td)c ~ 1. So, one gets, 

m (I /a) 

La ~ -p^- (32) 

Considering again the same parameters. Equation (|32l) gives k^a « 5 x 10^ for fif = 0.5, and 
k^a « 1.7 for yU/ = 0.95. One has to bear in mind that Equation (|3TI) is valid only for La «; 1, so 
we expect resonant absorption to be the dominant damping mechanism in the TT regime even for 
an almost neutral filament plasma. We will verify the analytical estimations of the present Section 
by means of numerical computations. 



2.4. Numerical Computations 

To numerically solve the full eigenval ue problem (Equations (|9l)-(fT4l)) we use the PDE2D 
code based on finite elements (jsewell 2005). We follow the same procedure as in Papers I and II. 



The integration of Equations (l9l)-(fT4l) is performed from the cylinder axis, r = 0, to the edge of the 
numerical domain, r = Vax> where all perturbations vanish since the evanescent condition is im- 
posed in the coronal medium, i.e., we restrict ourselves to trapped modes. The boundary conditions 
at r = are imposed by symmetry arguments. To obtain a good convergence of the solution and 
to avoid numerical errors, we need to locate the edge of the numerical domain far enough from the 
filament thread mean radius to satisfy the evanescent condition. We have considered Vax = 100a. 
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We use a nonuniform grid with a large density of grid points within the inhomogeneous transi- 
tional layer in order to correctly describe the small spatial scales that develop due to the Alfven 
resonance. 



3. RESULTS 



We focus our study on the ratio of the damping time to the period, td IP, which is the quantity 
that informs us about the efficiency of the kink mode damping. In the following sections, we 
compute Td//' as a fun ction of the dimensionles s longitudinal wave nurnber, kM . According to the 
observed wavelengths (|01iver & Ballesteiil2002h and thread widths (|Linl 120041) . the relevant range 
of k^a of filament oscillations corresponds to 10"^ < k,a < 10~^ So, the results within this range 
of k^a will deserve special attention. 



3.1. Filament Thread without Transitional Layer 

First, we start with the case of a homogeneous filament thread without transitional layer, i.e., 
l/a = 0. This is the configuration studied in Paper II with the addition of Hall's term in the 
induction equation. Hence, we briefly discuss these results here and refer the reader to Paper II 
for a more complete description. Figure displays t^/P versus k^a for different ionization 
degrees. In agreement with Paper II, we obtain that tu/P has a maximum which corresponds to 
the transition between the ohmic-dominated regime, which is almost independent of the ionization 
degree, to the region where Cowling's diffusion is more relevant and the ionization degree has 
a significant influence. The position of this transitional wavenumber is well approximated by 
Equation (fTSl) . As expected, the approximate solution obtained by solving Equation (|2TI) for ju = 
0.8 (symbols in Figure Wici)) agrees well with the numerical solution in the range of ha where 
Cowling's diffusion dominates, while it significantly diverges from the numerical solution in the 
region where ohmic diffusion is relevant. Within the range of typically reported wavelengths of 
filament oscillations (shaded region), t^/P is between 1 and 2 orders of magnitude larger than 
the observed values, which implies that neither ohmic nor Cowling's diffusion provide realistic 
damping times compatible with those observed. 

On the other hand, the kink mode exists as a propagating wave for k,a in the range k'l'a < 
ha < t^'^a, where h' and fc^"^ are the critical wavenumbers described by Equation (38) of Paper II. 
With the notation adopted in the present work, these critical wavenumbers can be expressed in 
dimensionless form as follows, 

* ? ^^^^y" ' - — ' (33) 

2 fjc 
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where kj^ is the perpendicular wavenumber to the magnetic field lines and contains the effect of the 
model geometry (see details in Paper II). For our present cylindrical filament thread, this geometry- 
related factor can be written as {k_i_df « {k,.df +{k^a^ , where kr and ky, are the radial and azimuthal 
wavenumbers, respectively. We approximate kj by its expression in the ideal case, 



CO] 



-kh 



A/ 



(34) 



which for ha 1 and pf simplifies to k^ ~ k^. On the other hand, the azimuthal wavenum- 
ber can be expressed as k'^ ~ rn^/a^. In the long-wavelength case, the azimuthal wavenumber 
dominates over the radial wavenumber, so we get {kj_a)^ ~ m^. Next, we plot in Figure |2tZ?) the 
ratio of the damping time to the period as a function of ha for /i/ = 0.8 considering three different 
values of the filament thread radius. We see that the smaller the thread radius, the smaller td/P 
and so the more attenuated the kink wave. In addition, the critical wavenumbers are shifted when 
a changes, the range of k^a for which the kink wave propagates being wider for thick threads than 
for thin threads. However, onemust bear in mind that the observed width of filament threads is in 
the range 0."2 - 0."6 (|Linl 120041) . and therefore a ranges from 75 km to 375 km, approximately. 
The critical wavenumbers are far from the relevant values of k^a for the observed thread widths, 
and so they should not affect the kink wave propagation in realistic filament threads. 




1000.00 








100.00 
















10.00 






'■■\\ 

'■\\ 

\\\ 


1.00 


- //.■ ' 
/ /.• 




\\ \ 
">\ \ 
\\ \ 1 




: / /•' 




'■A \ \ 




: / // 






\ 


0.10 


-i I : 
;) I ■ 






\ 

1 ^ 
1 : 


0.01 








i 


10'^ 10'* 


10''^ 


10° 10' 



k,a 



Fig. 2. — Ratio of the damping time to the period of the kink mode as a function of k,a correspond- 
ing to a thread without transitional layer, i.e., //a = 0. (a) Results for a = 100 km considering 
different ionization degrees: jif = 0.5 (dotted line), fif = 0.6 (dashed line), fif = 0.8 (solid line), 
and fif = 0.95 (dash-dotted line). Symbols are the approximate solution given by solving Equa- 
tion (EB for fif = 0.8. (b) Results for /// = 0.8 considering different thread widths: a = 100 km 
(solid line), a = 50 km (dotted line), and a = 200 km (dashed line). The shaded zone corresponds 
to the range of typically observed wavelengths of prominence oscillations. 



3.2. Filament Thread with an Inhomogeneous Transitional Layer 



Hereafter, we include an inhomogeneous transitional layer in the model. This causes the 
Alfven continuum to be present and so the kink mode is now resonantly coupled to continuum 
modes. In Figure [3ta) we assess the effect of the transitional layer on the ratio Tx^IP. For a fixed 
flf, we compare the results corresponding to several values of the transitional layer width. Some 
relevant differences appear with respect to the case Ija = 0. First, we see that t^/P is dramatically 
reduced for intermediate values of k,a including the region of typically observed wavelengths. 
In this region, the ratio t^/P becomes smaller as l/a is increased, a behavior consistent with 
damping by resonant absorption. This result is confirmed by means of Figures |4] and |5l which 
display the kink mode eigenfunctions close to the resonance point given by Equation (l23l) . We see 
that, with the exception of B,, which is directly proportional to the total pressure perturbation, the 
other perturbations show small spatial-scale oscillations of large amplitude around the resonance 
point, i.e., within the so-called dissipative layer. The azimuthal components of both the velocity 
and magnetic field perturbations have the largest amplitudes because of the coupling to torsional 
Alfven continuum rn odes. The thickness of the dissipative layer, 6a, is given by Equation (51) of 



By considering again that in our model the variations of the local Alfven frequency are only due 
to the variation of the equilibrium density, and approximating the resonant point by ~ a, we 
rewrite Equation (|35l) as follows. 



We see that 6a ~ l^^^, so the larger the transitional region width, the larger the thickness of the 
dissipative layer. This is consistent with Figures |4] and [5l in which the small spatial-scale oscilla- 
tions due to the resonance spread over a wider region for l/a = 0.5 (Figure |5]) than for l/a = 0.2 
(Figure SI). 

Turning back again to Figure Oa), we see that the ratio td/P for large k^a is independent of 
the transitional layer width and coincides with the solution in the absence of transitional layer. The 
cause of this behavior is that perturbations are essentially confined within the homogeneous part 
of the thread for large ha, and therefore the kink mode is mainly governed by the internal plasma 
conditions. On the other hand, the solution for small ha is completely different when l/a ^ 0. 
We note that the inclusion of the inhomogeneous transitional layer (i.e., for l/a 4^ 0) removes the 
smaller critical wavenumber, h^ , and consequently the kink mode exists for very small values of 



Sakuraietal.l(ll991ah . 
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(36) 
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Fig. 3. — Ratio of the damping time to the period of the kink mode as a function of k^a correspond- 
ing to a thread with an inhomogeneous transitional layer, {a) Results for = 0.8 considering 
different transitional layer widths: l/a = (dotted line), Ija = 0.1 (dashed line), l/a = 0.2 (solid 
line), and l/a = 0.4 (dash-dotted line). Symbols are the solution in the TB approximation given 
by solving Equation (l25l) for //a = 0.2. {b) Results for l/a = 0.2 considering different ionization 
degrees: fif = 0.5 (dotted line), fif = 0.6 (dashed line), p.f = 0.8 (solid line), and fif = 0.95 
(dash-dotted line). In both panels we have considered a = 100 km. 



(a) 



0-04 
0.02 
0.00 
-0.02 
-0.04 

-o.oe ^ 

-0.08 
-0.10 _ 



((>) 





1.0 




0.10 








-0.0000 










o.oa 




0.5 




-0.0001 
-0.0002 






0.5 




0.06 
0.04 










s 0.0 




°3 0.02 
0.00 




0.0 




-0.0003 






-0.5 


















-0.02 




-0.5 




-0.0004 






-10 




-0.04 








-0.0005 





1.0 
r/a 



1.0 
r/a 



Fig. 4. — Eigenfunctions of perturbations (a) v^, (b) v^, (c) B^, (d) B^, and (e) B^ of the kink 
mode around the resonant point (r/a « 1) for l/a = 0.2, k^a = 10~^, a = 100 km, and /// = 0.8. 
The solid line denotes the real part and the dotted line is the imaginary part of the corresponding 
eigenfunction. Arbitrary units have been used. 
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Fig. 5. — Same as Figure|4]but for l/a = 0.5. 



-17- 



We study in Figure [3tZ?) the dependence of tdIP on the ionization degree for a fixed 1 1 a, and it 
turns out that the ionization degree is only relevant for large k^a, where the ratio tq//' significantly 
depends of /iy . Figure |6] allows us to shed light on this result. There, we assess the ranges of 
k^a where Cowling's and Hall's diffusion dominate. As expected, we find that Hall's diff"usion is 
irrelevant in the whole studied range of k^a, while Cowling's diff"usion, caused by the presence 
of neutrals, dominates the damping for large k^a. In the whole range of relevant wavelengths, 
resonant absorption is the most eflicient damping mechanism, and the damping time is independent 
of the ionization degree as was analytically predicted in Section 12.3.11 On the contrary, we find 
that ohmic diffusion dominates for very small k-a. It is important noting that although the kink 
mode is still resonantly coupled to Alfven continuum modes, the damping time related to Ohm's 
dissipation becomes smaller than that due to resonant absorption, meaning that the kink wave 
is mainly damped by ohmic diffusion for very small k^a. Since rj is almost independent on the 
ionization degree, the ratio Tx^/P in the region of small k^a slightly depends on the value of /iy . 

Finally, we compare the full numerical solution with that obtained by solving the analytical 
dispersion relation in the TB approximation (Equation (|25l) ). For the case l/a = 0.2, and /i/ = 0.8, 
we plot by means of symbols in Figure Sa) the result obtained from Equation (|25] ). We can 
see a very good agreement between the approximation and the numerical result (solid line) for 
k^a > 10~^, while both solutions do not agree for ha < 10""^, which corresponds to the range of 
k^a for which ohmic diffusion dominates. Equation (|25l) was derived by taking into account the 
effect of resonant absorption and Cowling's diffusion, but the influence of ohmic diffusion on the 
damping is not included. For this reason, the approximate solution does not correctly describe 
the kink mode damping for very small k^a, which corresponds to extremely long and unrealistic 
wavelengths, while it successfully agrees with the full numerical solution for realistic wavelengths 
and larger values of k^a. 



4. CONCLUSION 

In this Paper, we have studied the combined effect of resonant absorption and partial ionization 
on the kink mode damping in a filament thread. Focusing on the results within the observationally 
relevant range of k^a, we have found that resonant absorption entir ely dominates the k ink mode 



damping, with tj)/P < 10 in agreement with previous investigations (lArregui et al.ll2008L Paper I). 
The obtained damping times are therefore consistent with those reported by the observations. The 
analytical value of Td /P in the TB and TT limits (Equation (|29l )) is a very good approximation to 
the numerical solution. None of the other dissipative effects studied here (i.e., ohmic, ambipolar, 
and Hall's diffusions) is of special relevance in the observed range of wavelengths, and the plasma 
partial ionization does not affect the mechanism of resonant absorption. The present results rein- 



10'^ 10" 10'^ 10° 10^ 

Fig. 6. — Ratio of the damping time to the period of the kink mode as a function of k,a correspond- 
ing to a thread with a = 100 km and l/a = 0.2. The different linestyles represent the results for: 
partially ionized thread with jj./ = 0.8 considering all the terms in the induction equation (solid 
line), partially ionized thread with fif = 0.8 neglecting Hall's term (symbols), and fully ionized 
thread (dotted line). 

force resonant absorption as the best candidate for the damping of filament thread oscillations. 

Some interesting issues not included here might be broached by future investigations. Some 
of them are, for example, to consider the longitudinal plasma inhomogeneity within the filament 
thread, to study the spatial damping of propagating waves, and to investigate the excitation and 
time-dependent evolution of impulsively generated oscillations. 
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